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ARRANGEMENT OF BODIES OF REVOLUTION IN SUPERSONIC 

FLOW TO REDUCE WAVE DRAG 
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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

ARRANGEMENT OF BODIES OF REVOLUTION IN SUPERSONIC 
FLOW TO REDUCE WAVE DRAG 
By Morris D. Friedman 


Summary 


The wave drag of a combination of slender bodies of revolution at 
zero angle of attack is studied with a view to determining the arrange- 
ments for which the total drag Is a minimum. Linearized theory is used 
to calculate the pressure distribution in the field surrounding the 
bodies. The interference drag coefficient is computed for different 
arrangements . 

The special cases of two bodies and of a three-body combination with 
b ila teral symmetry are considered. " The bodies treated are of the form 
determined by Sears and Haack as having minimum wave drag for prescribed 
volume and length. They also have equal fineness ratios. Numerical 
calculations of the drag coefficient of interference are carried out and 
curves are drawn which show the relative positions at which minimum drag 
occurs . 

A three-body configuration Is found for which the total wave drag 
Is about 35 percent less than the sum of the individual wave drags of 
the three bodies. 


INTRODUCTION 


The drag of a body of revolution in supersonic flow has been con- 
sidered by von Karman, Haack (references 1 and 2), and others. Haack 
and Sears (references 2 and 3) have determined theoretical body shapes 
for which the wave or pressure drag Is a minimum. Such bodies have 
important present-day flight applications since this wave drag, as an 
additional form of drag at supersonic speeds, limits the performance of 
modem aircraft . 
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drag is of interest. Different arrangements of the bodies of revolution 
at zero angle of attack are investigated, and that combination which has 
the smallest combined wave drag is determined. A practical combination 
would consist of a large body or fuselage and/or two smaller bodies 
which could represent either wing nacelles or tip tanks. In particular, 
the drag coefficient for the case of two slender bodies of equal fine- 
ness ratio (ratio of total length to maximum diameter) and unequal 
lengths is considered. (See fig. 1.) The procedure to be outlined, 
however, permits generalization, for the body shape chosen, to any num- 
ber of bodies of any relative size. 


SYMBOLS 


a speed of sound 

b abscissa of center of body of revolution 

c constant, the value of which determines the fineness ratio 

W 3 ) 

Cpj coefficient of interference drag 

Q]>t drag coefficient of the configuration based on the total fron- 

tal area of the configuration 

Cp pressure coefficient 

Dj drag change due to interaction of body pressure fields 

E complete elliptic integral of the second kind of modulus k 



F 0 f K 

K complete elliptic integral of the first kind of modulus k 

k modulus of the complete elliptic integrals 

L path of integration 

m half-length of body 

M Mach number 
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% free-stream dynamic pressure 

r cylindrical, coordinate (V y 2 + z 2 ) 

r 0 particular value of r 

R radius of a section 

RmflT radius of maximum section 

S cross-sectional area of body or bodies 

V free-stream velocity 

x,y,z rectangular coordinates 

0 Vm 2 - 1 

£ argument of the elliptic integrals of the third kind 

■Ao term which occurs in the circular case of the elliptic integral 

of the third kind 

£ dummy variable of integration 

<p perturbation velocity potential 

Subscripts 

parent body 
satellite body 
regions of integration 
combination as a whole 

GENERAL CONS LDERAPIOITS 

The type of configuration studied is illustrated in figure 1. It 
consists of a large and a small body each of the form having minimum 
wave drag for given volume and length and each of the same fineness ratio. 
The large body is situated on the axis with Its center at the origin. 

The smaller body has its axis parallel to the axis of the large body and 


1 

2 

I, II, III 
t 



4 


NACA RM A51I20 


may be. within the flow field of the large body. The equation for the 
shape of such bodies, which may be supposed to be created by distributing 
sources and sinks along an axis, is given by Haack (reference 2) as 

R 2 = c [m 2 - (x-b) 2 ] 3/a . (1) 


with b = 0 for the parent and b = b for the satellite body. (For 
b = 0, this formula is the expression for a thin body of revolution of 
length 2m with center at the origin and fineness ratio 2m/2R max . ) 

To calculate the interference drag of such a configuration, two 
possibilities must be considered: 

1. Only one of the bodies is within the disturbed flow field of 
the other . 

2. Each body is in the disturbed flow field from the other. 

When one of the bodies lies within the disturbed flow field of the 
other, the effects of the flow field in which this body is located must 
be considered. In other words. If all or part of either body is behind 
the Mach wave from the nose of the other body the streamlines will be 
distorted and a pressure will be exerted by the flow field of one body 
on the other. 

The potential field which results from the interaction of the flow- 
fields consists of the sum of the individual potentials and an Interfer- 
ence potential. In the cases when there are multiple reflections 
between the bodies a series of interference potentials may occur. Since 
the interference potentials are usually of higher order of smallness, 
they will be assumed negligible. Of course, at very high Mach numbers, 
or when the bodies are relatively close to each other, the effect of 
this interference may not be negligible. 

The interference drag coefficient due to the location of a body in 
the flow field of another can be evaluated by integrating the product of 
the additional, disturbed, pressure at a point and the slope of the body 
surface at that point. In a similar manner, the interference drag coef- 
ficient for the case when each body lies within the disturbed flow field 
of the other can be calculated. 

As a preliminary step, therefore, it is necessary to find an expres- 
sion for the pressure around each body for the whole region behind the 
Mach wave from the nose of the body. 
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METHOD OF CALCULATION 

Tne Pressure Field Surrounding A Single Body 


Under the assumptions of linearized theory, the shape of a slender 
body of revolution is described by the distribution along an axis of 
sources and sinks which satisfy the potential equation and the boundary 
conditions of uniform flow at infinity. Undter these assumptions the 
source strength is given by 


2sf’(x) = V- 


d% 

d3E2 


where S is the expression for the area of a section. 

The pressure coefficient, in this theory, is found from the relation 


„ 2 5 ? 

°T "fS 

where, as given in reference 2, 5cp/Sx is given by 

f , (6?ae 


S5 _ _ /*** 

3* 


'nose 


*/ (x-S) 2 - p^r 2 


( 2 ) 


(3) 


Here 


f(t) - Ja T p-“-- g Cl ^ 


Vn* - (t-b)= J 

is the expression for the source strength in terms of the following 
coordinates of integration: 


B 2 = M 2 - 1 
r 2 = y 2 + z 2 

When f(i) is replaced in equation (3) by its equivalent, and the 
substitution Is made for d <T/ox in equation (2), the pressure coeffi- 
cient is found to be determined by the integral 




[m 2 - 2(|-b) a ]di 


L J [m 2 - (t-b) 2 ][(x-5) 2 - B 2 r 2 l 


(h) 



6 


ID 


HACA EM A51I20 


This integral is elliptic, and the limits of integration and the method 
of solution depend on the following txrree regions (see fig. 2): 

Region X is hounded by the Mach aftercone from the nose and 
the Mach forecone from the tail. The limits of integra- 
tion are b-m and x-Br. 

Region II is bounded by the Mach forecone from the tail and 
the Mach aftercane from the tail. The limits of integration 
again are b-m and x-|3r . 

Region III is bounded by the Mach aftercone from the tail and 
infinity. The limits of integration . are b-m and b+m. 

With these limits the elliptic integrals for the different regions 
are complete and the solutions are: 1 

Region I 

C p = 


(m+3r) a - (x-b)‘ 


Em 2 + 23r(nH-b-x) - 2(x-b) 2 ] F 0 (ki) - 


[(m+3r) 2 - (x-b) 2 ] E 0 (ki) + 

2(x-b) J (m+pr) 2 - (x-b) 2 A Q (kj, Sj) J>- (5) 

where 



transforms equation (4) into normal elliptic form. 




elliptic integrals compiled by Mir. Paul F. Byrd of the 
Ames Aeronautical Laboratory. 
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Region II 


C P = — ~ ^ 7C — Tm(m+2pr+2b-2x)F 0 (kii)-lfm3rE 0 (k II )+4(x-b) «/ mgr £n)|- 

2i/mgr L (gj 


where 


. / (*-bf - (^r)* ; 

y l®Br y x-b+| 


2m 

•b+'gr+m 


and the change of variable 


-1 / 2 

u = an 1 / 

V (x-b- 


2gr(£-b-Hn) 


-gr+m)(x+gr-i) 
transforms equation (4) into normal elliptic form. 

Region III 

- 7 . ===== T (x-b-prJ^oCinx) + [ (x-b ) 2 - (m-gr) 2 ] 

V(x-b ) 2 - (m-gr ) 2 L 

^(^IIl)! ” 2(x-b ) J (x-b) 2 - (m-gr ) 2 Ao(k TTT , C TT j)j- 


Cp = 


where 


^III 


/(x-b 


ij-Brm 


(x-b) 2 - (m-gr) 2 
and the change of variable 


f . /x+gr-m-b 

> tEII = sin -l ~ 7 - 

V x+gr-m-b 


( 7 ) 


u = sn- 1 /U+gr-m-b)(S-b-m) 
v 2m(x+gr-§ ) 


transforms equation (4) into normal elliptic form. 
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The functions Fo, Eo> and Ao are tabulated in reference 4, where 
they are also defined as: 

FoOO = ^K(k) 

E 0 (k) = -f-E(k) 

Ao(k^) = E 0 (k)F(k',5) + F 0 (k)E(k*,5) - F 0 (k)F(k»,S) 


where 

K(k) complete elliptic integral of the first kind of modulus k 


E(k) 


complete elliptic integral of the second kind of modulus k 


F(k*,0 

E(k*,£) 


} 


incomplete ellipti c int egrals of the first and second kind of 
modulus k* (= 1-k 2 ) , and argument £, respectively 


Thus, according to linearized theory, equations (5), (6), end (7) 
completely deter min e the pressure field in the region behind the Mach 
wave from the nose of a body of revolution of the prescribed shape. An 
isometric sketch of the pressure coefficient at fixed values of the 
radial distance is shown in figure 3« 


Interference Drag Coefficient 


An inspection of figure 3 shows that the gradient of the curve of 
the pressure coefficient in the stream direction changes from negative 
to positive behind the center of the body. This region of positive 
pressure gradients is a zone of favorable buoyancy in the pressure field 
around the body. Therefore, a small body placed anywhere within this 
zone should have a negative interference drag (i.e., a thrust) due to 
this pressure field which should cause a decrease in the combined drag. 

If the pressure coefficient surrounding a body at a given lateral 
distance r = r 0 be called Cp^, the interference drag coefficient may 
then be calculated from the expression 


_ 1 r tail _ _ RdR , 

Dl *Rmax 2 Jnose Pt ^ 


-'nose 


(8) 
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where Rmax is the radius, of the frontal area of the reference body. 
Substitution of the respective values for the quantities in equation (8) 
yields 

QDj = '~^ av £ C Pt{^ _3,tC ( x ” b ) - (x-b) a j- to (9) 

Actually, the integral in equation (8) or (9) is made up of one or more 
integrals depending on the number of regions in which the small body is 
located. The limits of integration depend on the parameters b, m, r, 
and 3 which determine the regions of integration. This integral must 
be computed numerically. 


For a clearer exposition of the mathematical computations in the 
cases where interference exists, it is necessary to distinguish between 
the bodies. In the present case, if the central body situated on the 
x axis be denoted by the subscript 1 and the other body by the sub- 
script 2 , then, depending on which body is being acted upon, the inter- 
ference drag coefficient is 


pb z +m z s- 

CDi = — 2 / Cp t J-3Jtc 2 (x-b 2 ) [m 2 2 - (x-b 2 ) 2 ] 1/2 t-to (10) 

tfRmax '-'bo-m.a 1 J 


or 


itRmax 1 -'b 2 -m 2 


obi+mi f 

C »z = °Pt a {-3*=xC-»x) - (x-bx) 2 ] 1 ^ (11) 


Equations (10) and (ll) are, to a first approximation, the drag 
coefficients due to the effect of body 1 on 2 and body 2 on 1, respec- 
tively. The integration extends from the nose to the tail of the body 
which is in the disturbed flow field. Similarly, for the calculation, 
to a first approximation, of the drag coefficients due to the inter- 
action of the bodies, formulas (l0) and (ll) are both used. 

The preceding equations permit the calculation of the interference 
drag coefficient for any arrangement of bodies. In particular, equa- 
tions (10) and (ll) apply to the special- case when two identical bodies 
are symmetrically placed with respect to a central body so that each 
body may be in the flow field of the others. The formulas may be applied 
to any configuration, whatever the interference pattern. 
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Numerical Calculations 


As an exploratory investigation, the simple case of a small body • 
having one-half the length of the large body was chosen. This combina- 
tion is derived by assigning the following values to the parameters in 
equation (l): 


c m b 

Ci = 0.005 mi = 2 bi = 0 

c 2 = .01 m 2 = 1 b 2 = b 

The equations which describe the bodies are 

R x 2 = 0.005 (ii--x 2 ) S/2 

R 2 2 = .01 [ 1 - (x-b) 2 ] 3/2 

It is to be noted that each body has the fineness ratio 10 but that one 
is twice as long as the other. 

The pressure coefficient depends on 0 and r 0 . The dependence Is 
on 0r o , rather than on either parameter individually, and the calcula-. 
tions were carried out for 0r o = 0.5m 2 > l.Cm 2 , and 2.Cta 2 ; where,, as 
shown previously m 2 =1. A change in 0 and in r 0 which keeps 0r o 
invariant will not change the calculations . 


The drag coefficient, with respect to the frontal area, of the 
individual body is given by Haack (reference l) as 

A AQ-. 2 

■ The interference drag coefficient is computed from 


or in 


this case 


8 


formula ( 9 ). For each value of 0r o , therefore, the remaining parame- 
ter b, in this case the longitudinal distance between the centers of 
the bodies, is given values between -3m 2 and +1 hh 2 . 


The interference drag coefficients based on frontal area of the 

small body are sketched for different and JL in figure 4. Since 

m 2 m.2 


the reversibility of drag (reference 5) holds true, the curves would be 

"b b 

expected to be symmetrical about the line — = 0 (where ■ — is the 

I II 2 Blp 

nondimena ional longitudinal distance between the centers of the bodies). 
The symmetry in the figure is a measure of the accuracy of the numerical 
calculations. From the figure. It Is evident that minimum drag occurs 

when the small body Is close to the large one ( ^£2. = 0.5/ 

\ m 2 / 

Its center just forward of the tail wave from the large body. 


and with 
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It may be observed that, while the interference drag varies markedly 
with fore and aft position, for the range of j3r 0 considered the maxi- 
mum. favorable interference that can be obtained gener a lly decreases with 
increasing separation between the bodies. It is interesting to note 
that drag min imums occur whenever the center of the satellite body is 
Just forward of the tail Mach wave of the parent body. Also of interest 
is the fact that, in the case when both bodies had their centers on the 
vertical axis, the lowest drag occurred when the lateral separation 
parameter ( pro) was equal to the length of the small body and higher 
drags resulted as the bodies were brought closer together. 

In figure 5, there are plotted the values of the total, drag coeff i- . 
cient, based on the total frontal area of a three-body combination with 
bilateral symmetry, against the longitudinal distance between the centers 
of the parent and satellite bodies. Since interaction between all three 
bodies occurred only for the lowest value of gr 0 (= 0.5m 2 ) where it was 
found to be negligible, the interference drag coefficients previously 
calculated for the two-body configuration could be used directly to 
determine the total drag of the three-body configuration. As a conse- 
quence, the variations in drag coefficient are similar to those of 
figure 4. Again the lowest drag occurred at pro = 0. and. was 
approximately 35 percent less than the drag of the three bodies without 
interference. 

Because of the -unusual shapes of the curves in figure 4, it was 
considered advisable to investigate the drag interaction of a combina- 
tion of bodies of different shape. The slender pointed body derived by 
Jones and Margolis (reference 6) was selected for this purpose since for 
the. same fineness ratio the drag coefficient is comparable to that of 
the Sears -Haack body. The interference drag coefficients of a combina- 
tion of two such bodies with m 1 /m 2 = 2 were calculated for b = 0, and 
different values of gr 0 . Reasonable agreement was obtained with the 
results of figure 4 for pr Q = 2.Ctn 2 and £3r 0 =' 1.0m 2 . In the case when 
pr Q = O.^m. 2 , since the bodies are close together, there Is a discrepancy 
which may he due to the differences in body geometry. 


CONCLUDING^ REMARKS 


It Is found that the combined wave drag of a combination of bodies 
of revolution can be decreased if an arrangement is chosen which takes 
advantage of a favorable pressure zone which exists behind the center of 
each body. 

In the cases discussed, where the ratio of the lengths of the 
bodies was 2 to 1 but the fineness ratios were equal, numerical calcula- 
tions showed that the maximum favorable interference occurred when the 
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center of the small tody was just forward of the stern Mach wave from 
the large tody and the todies were close together. For the range of 
|3r 0 considered the magnitude of the favorable interference generally 
decreased with Increasing separation between the todies. 

In the case of a bilaterally symmetrical arrangement of three bodies 
with a lateral separation equal to one-quarter of the length of the small 
body the total wave drag was found to be 35 percent less than the com- 
bined wave drag of the three bodies. 


National Advisory Committee for Aeronautics 
Ames Aeronautical Laborato^*y 

Moffett Field, Calif., Sept. 20, 1951 
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Figure /.- Coordinate system. 
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Figure 5*- Total drag coefficient of three-body combination with bilateral symmetry. 
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